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Abstract 

Spin dynamics of a polarized spin system is studied when the latter 
is coupled with a resonant electric circuit and is under the action of 
an external pumping supporting a stationary nonequilibrium magne- 
tization. A complete classification of possible regimes of spin motion 
is given. In addition to seven regimes considered earlier, two other 
transient regimes are found and thoroughly described: One is an os- 
cillatory regime, when spins always move coherently but the degree 
of coherence fluctuates with time. Another is a pulsing regime, when 
spins reveal coherent motion during short pulses separated from each 
other by intervals of incoherent motion. These regimes are, in princi- 
ple, transient, although may be extremely long lasting; their duration 
may be several orders longer than the transverse relaxation time and 
twice longer than the longitudinal relaxation time. Both transient 
regimes end with a coherent quasistationary regime. 
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1 Introduction 



Nonequilibrium resonance phenomena in spin systems have their counter- 
parts in atomic systems. Recall, for instance, free induction, occurring sim- 
ilarly in both kinds of the systems, or spin echo being a direct analog of 
photon echo. The feasibility of a self-organized coherent process, called su- 
perradiance, has been theoretically predicted almost simultaneously for spin 
[1] and atomic [2] systems. The difference between these is that for real- 
izing such a coherent process in a spin system, the latter must be coupled 
with a resonator. Superradiance in atomic and molecular systems has been 
studied, both theoretically and experimentally, quite in detail, and has been 
expounded in a number of reviews and books, of which we cite only some 
recent Refs. [3-6]. 

In analogy with atomic superradiance, the process of collective coherent 
relaxation in spin systems has been called spin superradiance. This process 
was observed for electron [7-9] as well as for nuclear spins [10-12]. Accu- 
rate experiments observing purely self-organized superradiance from proton 
spins have been accomplished [13-16]. The peculiarities of spin superradi- 
ance were studied by means of computer simulation [17,18], being based on 
the microscopic Hamiltonian of a nuclear magnet, commonly accepted in the 
theory of nuclear magnetic resonance [19]. An analytic theory for this mi- 
croscopic model was developed in Refs. [20,21], where it was shown that, 
despite many similarities between spin and optic superradiance, there are 
also crucial differences between them. For instance, because of the principal 
role in triggering the radiation process, played by direct dipole interactions, 
the self-organized coherent relaxation in nuclear magnets is of non-Dicke 
type [21,22]. The behaviour of resonant spin systems under the action of 
an injected signal with particularly chosen delay times [23] and under the 
influence of parametric excitations [24-26] has also been considered. The im- 
portance of studying nonequilibrium coherent phenomena in spin systems is 
caused by the usage of these phenomena for many applications, for example, 
for spin masers [27-30], for the repolarization of scattering targets [16,31], 
and for a possible creation of sensitive particle detectors [32]. 

There is a problem that has not yet been properly studied for nuclear spin 
systems coupled with a resonator: What would be the behaviour of such a 
system under the influence of a stationary nonresonant pumping supporting 
a nonequilibrium magnetization? This kind of pumping could be realized by 
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means of dynamical nuclear polarization. Note that an equivalent question 
was posed for electron spin systems [33] and considered in the framework of 
the adiabatic approximation. However, the latter, as is well known [34], is 
valid only at the last stage of relaxation, when the system is already close to 
its stationary state. The adiabatic approximation cannot describe a transient 
process, as has been discussed in detail in Refs. [21,35]. Therefore, the 
authors of Ref. [33] considered only the asymptotic stationary regime for 
a spin system with a constant nonresonant pumping. Such a problem for 
atomic systems has been analysed earlier and it has been shown that in 
the presence of a constant external pumping atomic systems may exhibit 
only pulsed operation and cannot work in the quasistationary regime (see 
discussion in Refs. [36,37]). 

In the present paper, we consider a nuclear magnet coupled with a res- 
onant electric circuit and subject to the action of a constant nonresonant 
pumping supporting a nonequilibrium stationary magnetization. The con- 
sideration is based on the standard microscopic Hamiltonian [19] with dipole 
interactions between nuclear spins. We do not invoke the adiabatic approx- 
imation, but use the scale separation approach [20,21]. Therefore, we may 
analyse all possible transient regimes. 



2 Nuclear Magnet 

The standard Hamiltonian modelling a solid sample consisting of N nuclear 
spins can be written [19] in the form 
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H = ~T, H a-»T,B- Si, 

z i+j i=l 



(1) 



in which 



ij 



Si ■ Sj 3 [ Si ■ n>a ) I Sj ' n> 



is the dipole interaction energy; //, a nuclear magneton; Si— {Sf, Sf, 5?}, a 
spin operator; and 



fij — | f ij \i — 



_ r ij 



r ij— r i r j ■ 
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The total magnetic field 

B=H + H, H =H e z , H=He x , (2) 

consists of an external magnetic field Ho directed along the z axis and of a 
feedback field H of the coil of a resonator electric circuit with the coil axis 
being directed along the axis x. The resonator magnetic field H is formed 
by an electric current satisfying the Kirchhoff equation. 

Let us define the Larmor frequency u and the resonator frequency u, 
being, respectively, 

where Ho is assumed to be positive, L is the coil inductance, and C is the 
circuit capacity. Introduce also the resonator ringing damping. 

in which Q is the quality factor of a circuit and R, its resistance. Then the 
Kirchhoff equation can be written in the form 

dH 9 ft TT . ... A dM x 

where r\ is the coil filling factor; p, the density of spins; and 

a N 
i=i 

is the x component of the reduced magnetization. 

There are the following characteristic times: the spin-lattice relaxation 
time Ti; the spin-spin dephasing time T 2 ; the inhomogeneous broadening 
time T 2 * , due to local random fluctuations; and the resonator ringing time 
T 3 . The width corresponding to these times are 

= J_ = J_ * = — - — 

Tl rp ) ^2 rp 1 12 rp* ? 73 rp • 

J l J 2 J 2 J 3 

In the presence of pumping with the pumping velocity j p , the effective lon- 
gitudinal relaxation becomes 

7i =7i + 7 P - ( 6 ) 
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All these widths are usually small, as compared to the frequencies in (3), 
defining the set of small parameters 

^-<1, — <1, ^<1, — <1. (7) 

C^o ^0 ^0 ^ 

The resonator natural frequency is tuned to be close to the Larmor frequency, 
so that the detuning from the resonance is small, 



|A| 



< 1, A = u-cu . (8) 



The existence of small parameters in (7) and (8) justifies the use of the 
scale separation approach, whose all details have been thoroughly expounded 
in Ref. [21]. Employing this approach, we may derive the evolution equations 
for the transverse magnetization 



1 N 



«^E<^-tf?> (9) 



and the longitudinal magnetization 



iV 3=1 

The coupling between the spin sample and the resonant electric circuit is 
described [20,21] by the effective coupling parameter 

It is convenient to introduce the function 

w = v 2 z 2 , v= \u\. (12) 

\LUoJ 

After averaging over the time 2it/ojQ of fast oscillations and over random local 
fields, we obtain [21] the system of equations for function (12), 

dw , , . 

— = -2 l2 (l + gz)w, (13) 



and for the longitudinal magnetization (10), 



dz 
~dl 



g l2 w - 7* (z - C) ■ 



(14) 



The derivation of eqs. (13) and (14) has been carefully explained in Ref. 
[21]. The only difference, in the case we consider now, is that the spin- 
lattice relaxation constant 71 is replaced by the effective longitudinal width 
7* — 7i + 7 P , including the pumping velocity j p , and that the stationary 
magnetization parameter (, in the presence of pumping, becomes negative, 
C < 0. The evolution equations are complemented by the initial conditions 



In what follows we assume that the coupling parameter (11) is nonzero, since 
the case g — > would result in the trivial exponential relaxation of solutions 
to Eqs. (13) and (14). 

The spin-lattice relaxation parameter 7x is usually much less than the 
dephasing width 72, although for some materials they can be rather close to 
each other. For instance, in the case of z He at low temperature T ~ IK, 
for the characteristic times one has [38] the values 7\ = 100 — 300s and 
T 2 = 30 — 100s, so that 71/72 ~ 1/3. And as far as the effective relaxation 
parameter in the presence of pumping is the sum 7* = 71 +7 P , hence 7* > 71, 
as a result of which 7* can become comparable with 72 even if 71 <C 72. 

What is more significant is that, even, when 7^ is negligibly small as 
compared to 72, the term in (14) containing 7^, as will be shown in what 
follows, cannot be omitted if £ < 0, that is, if a pumping is present. This 
situation is drastically different from the case when pumping is absent [20- 
22]. In the latter case, if 7^ <C 72, then in the time interval < t < (7i) _1 
one can omit the last term in (13). This omition allows to solve the system 
of Eqs. (13) and (14) exactly, yielding 



w(0) = w , 



z(0) = z . 



(15) 




(16) 



and 




(17) 
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where the radiation width is 



7o = 72 



(18) 



the radiation time is r = 7 1 , and the delay time is 



Jo_ ln 7o ~ 72(1+^^0) 

2 70 + 72(1 + gz ) 



(19) 



Solutions (16) and (17), depending on initial conditions and the coupling 
parameter (11), describe seven qualitatively different regimes of spin relax- 
ations: free induction, collective induction, free relaxation, collective relax- 
ation, weak superradiance, pure superradiance, and triggered superradiance. 
This classification is valid when the pumping is absent, for which case all 
these regimes have been analyzed earlier [20-22]. 

In the presence of pumping, we may expect that solutions (16) and (17) 
correctly describe the beginning of the relaxation process for time t <C (7i) _1 - 
These solutions, when < 7^0 < 1) depict the first superradiant burst 
occurring at the time t — to, where 



However, these solutions do not give the overall picture for all times, even 
when 7^ <C 72- As the analysis of the following sections shows, the whole 
behaviour of solutions to Eqs. (12) and (13), in the presence of pumping, is 
essentially more complicated. 

3 Stability Analysis 

In order to understand, from the mathematical point of view, why the solu- 
tions to Eqs. (13) and (14), if one puts 7^ zero, can be drastically different 
from those when 7^ 7^ 0, even if 7* is negligibly small as compared to 72, 
one has to accomplish the stability analysis. For this purpose, we write Eqs. 
(13) and (14) in the form 




(20) 





(21) 
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in which 



-272(1 + gz)w, V 2 = g~f 2 w - 7*(z - C)- 



The equations Vi = V 2 = define stationary, or fixed, points. We have two 
such points, one is 

Z * = c, wi = 



and another is 



4 = — , w 2 = — Y~( 1 + gC) 

9 9 2 l2 



The Jacobian matrix 





- dVi 


dVi - 




dw 


dz 


J = 








av 2 


8V 2 




- dw 


dz J 



(22) 
(23) 

(24) 



corresponding to (21), takes the form 

-272(1 + 5^) -2^ 2 gw 



J = 



912 



-11 



(25) 



The eigenvalues of matrix (25) are 
1 



7* + 2 72 (1 + gz) ± (7* - 2 72 (1 + gz)Y - Sg 2 ^w . (26) 



1/2 



The values of (26) evaluated at the fixed points define the Lyapunov expo- 
nents. At the first fixed point, given by (22), we have 



K = ~ll Ar = -2 72 (1 + <7C)- 
And at the second fixed point, given by (23), we find 

^2 = ~\ \ll ± ^(ll) 2 + 87l*72(l+5C)} • 



(27) 



(28) 



The stability of the fixed points and, consequently, the stability of motion 
is characterized by the signs of the Lyapunov exponents [39]. Varying the 
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value of the pumping parameter £, we may get qualitatively different regimes 
of motion. These regimes are separated by the pumping thresholds 

C^-i C^-Ul + f). (29) 
9 9 \ 872/ 

When the pumping parameter satisfies the inequality 

C > G, (30) 

then 

Af < 0, A^ < 0, A^ > 0. 

Hence, the fixed point (22) is a stable node, while that (23) is a saddle point. 
In the case when 

C = G, (31) 

we have 

A+ = A+ < 0, \~ = \ 2 = 0. 

Both fixed points (22) and (23) merge together becoming neutrally stable. 
The pumping threshold £i corresponds to a bifurcation point. 
When the pumping parameter is in the region 

C2 < C < G, (32) 

then 

A+ < 0, Af > 0, A^ < 0, 

which means that the fixed points interchange their properties: now (22) is 
a saddle point and (23) becomes a stable node. 
For sufficiently strong pumping, when 

C < C2, (33) 

the fixed point (22), as earlier, continues to be a saddle point, since < 
and Af > 0. But for the fixed point (23) we get 

A£ = ~7r±<n, (34) 
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with 

« = \ \/l(7!) 2 + 871*72(1 + ^)1- (35) 

This shows that (23) transforms to a stable focus. 

Let us notice that under all conditions, if one puts 7* — > 0, then (27) and 
(28) yield A^ = and \f = 0. This tells us that both fixed points correspond 
to a structurally unstable system [39]. Structural instability means that 
the temporal behaviour of the system can be essentially disturbed under an 
arbitrary small change of the evolution equations. This is why the term of 
Eq. (14) containing 7*, in general, must not be omitted even if 7^ is many 
orders smaller than 7 2 . It may happen that for time t <C (7*) -1 , one can 
omit this term in some particular cases. For instance, this is the case of a 
spin maser without pumping [20,21], that is, with ( > 0. However, for a 
spin maser in the presence of pumping, when ( < 0, the situation can be 
drastically different. In the latter case, to describe the temporal behaviour 
of the system, one has to keep the term with 7*. 



4 Numerical Solution 

To analyse the behaviour of solutions to the system of equations (13) and 
(14), we solved this system numerically. The spin-resonator coupling param- 
eter (11) is taken to be g = 10. For a weak pumping, when ( > (1 = —0.1, 
as well as for an intermediate pumping, when ( satisfies (32), the behaviour 
of solutions is similar to that studied in Refs. [20-22,28-30], which is caused 
by the fact that the stationary point in all these cases is a stable node. The 
most interesting here is the case when the fixed point is a stable focus. Then 
qualitatively new types of solutions appear. Therefore, in what follows we 
concentrate our attention on the case of strong pumping corresponding to 
inequality (33). To this end, we take the pumping parameter ( = —0.5. A 
few typical phase portraits for a system with a fixed point being a stable 
focus are presented in Fig.l. For convenience, we introduce a notation 




Increasing the pumping velocity leads to the increase of 7, as a result of 
which the stationary value of vol a l so increases. Note that the initial value 
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of zq = z(0) does not influence much the whole picture. The phase portraits 
for z = —0.5 and z = +0.5 are very similar to each other. 

The following figures show the temporal behaviour of the functions w(t) 
and z(t) for various initial conditions and pumping velocities. When there 
are no external pumping fields, except that realizing a stationary dynamical 
polarization (, then we should put z(0) = (. However, the initial polarization 
Zq = z(0) can be made different from ( by using additional short external 
pulses. Keeping this possibility in mind, we consider the cases with z(0) not 
always coinciding with (. 

In Figs. 2 to 6, we see the oscillatory regime of motion. Everywhere, 
if it is not stated otherwise, we take g — 10 and ( = —0.5. It is only in 
Fig. 7 where the pumping parameter is varyed. Fig. 8 demonstrates how an 
oscillatory regime of motion changes to a pulsing one with changing initial 
conditions. The pulsing regime of motion is presented in Figs. 9 to 12. 

5 Discussion 

We have considered the dynamics of spin maser under a stationary pumping 
supporting a constant nonequilibrium magnetization in a system of nuclear 
spins. Such a pumping can be accomplished by means of stationary dynam- 
ical polarization of nuclei. The regimes of oscillatory and pulsing motion are 
found. The distinction between these two regimes is, of course, somewhat 
conditional, depending on the level of coherence existing during the time 
intervals separating the neighbouring coherent bursts. These bursts can be 
directly observed by measuring the current power P that is proportional to 
v 2 . As far as 7^ C w , we have v 2 ~ w. Whence, P ~ w. In this way, the 
function w(t) is proportional to the current power and, thus, is an observable 
quantity. This function is also connected with the intensity of magnetodipole 
radiation, I(t), and the coherence coefficient C co h(t), as is discussed in Refs. 
[17,18]. Therefore, the coherent bursts occurring in the system of nuclear 
spins can be named superradiant. 

Superradiant regimes appearing in a spin maser under the action of the 
nonresonant pumping, supporting a stationary level of a nonequilibrium mag- 
netization (, are quite different form the regimes developing under the influ- 
ence of a resonant pumping realized by means of an alternating external fields. 
For comparison, we adduce in Figs. 13 to 15 the behaviour of the coherence 
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coefficient C co h, radiation intensity /, and of the polarization p z = —z(t) for 
a nuclear magnet pumped with a resonant alternating field [17,18] oscillating 
as h cos cut. 

Dynamics of a spin maser in the presence of a nonresonant pumping 
supporting a pumping polarization ( < ( 2 resembles that of pulsing lasers 
[36,37]. The difference is that pulsing lasers cannot operate in a stationary 
regime, while a pumped spin maser, after its oscillatory or pulsing stage, 
tends to a stationary regime with a current power proportional to w 2 in (23). 
The level of coherence in this stationary regime, as compared to that of the 
first superradiant burst, defined in (20), is described by the ratio 

W*2 _ 7i|l+flCI 

w(t ) 72(1 + gzo) 2 + g 2 wo' 

The latter, for w = 0, z ~ ( ~ 1, and g 3> 1, gives 

w(t ) 7 2 ff" 

Therefore, to reach the intensity of the first superradiant burst, the pumping 
velocity is to be very high, so that 7* ~ 7 2 g which looks as practically 
unattainable. 

It is worth paying some attention to terminology. The temporal behaviour 
of solutions, as is seen from above figures, is not periodic, since the time in- 
tervals between pulses as well as their amplitudes change with time. This 
time dependence cannot be called quasiperiodic (in mathematical sense). It 
is not chaotic too, although for some parameters and in a limited time inter- 
val it may look as pseudochaotic, slightly reminding quantum pseudochaos 
[40]. Therefore, the most appropriate adjectives characterizing the type of 
solutions we found could be, probably, oscillating and pulsing. Or we could 
describe all of them by one word, e.g., pulsing, implying that this incapsu- 
lates all admissible variants of solutions consisting of a number of pulses. 
Such a pulsing operation can be employed in spin masers. 
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Figure Captions 



Fig.l. The phase portraits for g — 10, ( = —0.5, w = 0,001, and 
for the time interval < t < 1007^" 1 for different pumping velocities and 
initial polarizations: (a) 7 = 0.1, Zo = 0.5; (b) 7 = 1, z = —0.5; (c) 
7 = 1, z = 0.5. 

Fig. 2. The temporal behaviour of solutions for 7 = 0.001, with initial 
conditions w = 10~ 6 and z = —0.1: (a) w(t); (b) z(t). 

Fig. 3. Solutions to evolution equations for 7 = 0.001, with initial condi- 
tions w = 10 -6 and z = —0.1: (a) w(t); (b) z(t). 

Fig. 4. Solutions of equations as functions of time for 7 = 1 and initial 
conditions w = 0.001 and z = —0.5: (a) w(t); (b) z(t). 

Fig. 5. Spin dynamics for initial conditions wo = 0.1 and z$ = —0.25, with 
different parameters 7, where the solid line is for 7 = 1 while the dashed line 
is for 7 = 0.5: (a) w(t); (b) z(t). 

Fig. 6. Spin dymanics with initial conditions w = 0.5, z = 0.5 for 7 = 1 
(solid line) and 7 = 0.5 (dashed line): (a) w(t); (b) z(t). 

Fig. 7. Spin dynamics for 7 = 1 with initial conditions wo = 0.5, z = 0.5 
and a varying pumping parameter ( = —0.5 (solid line), ( = —0.3 (dashed 
line): (a) w(t); (b) z(t). 

Fig. 8. Transformation of an oscillatory regime of motion for w(t) to a 
pulsing one, under 7 = 0.01, when changing initial conditions: (a) w = 
0.001, z = -0.1; (b) w = 0.01, z = -0.1. 

Fig. 9. Pulsing regime of motion for 7 = 0.1, with initial conditions 
w = 0.01 and z = -0.5: (a) w(t); (b) z{t). 

Fig. 10. Pulsing regime of motion for 7 = 0.01, with initial conditions 
Wo = 0.01 and z = 0.5: (a) w(t); (b) z{t). 
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Fig.ll. Change in the behaviour of the function w(t), with the same 
initial conditions wo = 0.001, z = —0.5, under the variation of the effective 
longitudinal relaxation: (a) 7 = 0.1; (b) 7 = 0.001. 

Fig. 12. Dynamics of the function w(t) in the pulsing regime with dif- 
ferent parameters: (a) 7 = 0.1, Wo = 10~ 6 , z = —0.5; (b) 7 = 0.01, wo = 
0.1, z = -0.1. 

Fig. 13. The coherence coefficient C co h, intensity of radiation / in arbi- 
trary units, and polarization p z = —z versus time for g = 0, ujq = 2OO72, and 
uj = 2OO72. The influence of different intial conditions is analyzed: zq = 0.475 
(solid line); z = —0.375 (dashed line); z Q = —0.475 (solid line with crosses). 

Fig. 14. The same functions as in Fig. 13 for a spin system with switched 
off dipole interactions in the case of g = 0, ujq = 2OO72, and z = —0.475. The 
external alternating field is not in an exact resonance, with the frequencies 
uj = IOO72 (solid line) and ui = 19572 (dashed line). 

Fig. 15. The coherence coefficient C co h and radiation intensity / versus 
time for g = 0, c^o = 2OO72, uj = 20572, and Zq = —0.475. Different number 
of spins is considered: N = 343 (solid line), N = 125 (dashed line), and 
N = 27 (solid line with squares). 
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